Landau-Zener Transitions in an Adiabatic Quantum Computer 
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We report an experimental measurement of Landau-Zener transitions on an individual flux qubit 
within a multi-qubit superconducting chip designed for adiabatic quantum computation. The 
method used isolates a single qubit, tunes its tunneling amplitude A into the limit where A is 
much less than both the temperature T and the decoherence-induced energy level broadening, and 
forces it to undergo a Landau-Zener transition. We find that the behavior of the qubit agrees to a 
high degree of accuracy with theoretical predictions for Landau-Zener transition probabilities for a 
double-well quantum system coupled to 1/f magnetic fiux noise. 
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Adiabatic quantum computation (AQC) is a quantum 
mechanical method for solving hard computational prob- 
lems [1[. In AQC one encodes a computational problem 
in a suitable physical system. If one can somehow put the 
system in its ground, or lowest energy state, the structure 
of that ground state then reveals the answer to the prob- 
lem. To find that ground state using AQC, one starts by 
engineering a simple Hamiltonian or energy functional for 
the system, and by placing the system in the ground state 
of this simple Hamiltonian. Then one gradually deforms 
the Hamiltonian of the system from the simple form into 
a complex Hamiltonian whose ground state encodes the 
answer to the problem. If this deformation is sufficiently 
gradual, then the transformation of the state of the sys- 
tem is adiabatic, and the system remains in its ground 
state throughout the deformation. AQC is known to be 
a universal model of quantum computation [2| . 

A key question of AQC is whether adiabaticity can be 
maintained throughout the computation. For an isolated 
system, Landau-Zener (LZ) transitions [3|, Ij] may take 
the system out of its ground state, unless the time over 
which one transverses the point of the minimum energy 
gap during the course of a computation is longer than 
the instantaneous coupling between ground and excited 
states divided by the minimum gap squared. Just which 
hard problems can be encoded in such a way that adi- 
abaticity can be maintained over reasonable times is an 
open question. Moreover, for some optimization prob- 
lems, an excited final state with sufficiently low energy 
could provide an acceptable solution and therefore the 
adiabaticity condition may be relaxed. 

Equally important question is whether interactions be- 
tween the computer and its environment can spoil the 
computation. Unlike in the gate model of quantum com- 
putation i5|, the effects of an environment on AQC are 
less understood. While it is now clear that AQC has fun- 
damental advantages over the gate model in regards to 
robustness against decoherence [a [Zl [sj, |9| , there does not 
yet exist an equivalent of the threshold theorem [5| that 
describes under what conditions AQC coupled to an en- 
vironment will succeed. Nonetheless, a clear understand- 



ing of how LZ transitions are affected by an environment 
represents an important step forward. 

For large-scale hard problems, it is inevitable that any 
implementation of AQC will encounter minimum gaps 
that are smaller than temperature T and decoherence 
rate. In order to understand the effects of environment 
in this limit, we attempt to first understand in detail how 
environment affects LZ transitions in individual qubits 
within an adiabatic quantum computer. We isolate a 
single qubit in a 28 qubit superconducting chip from its 
surrounding qubits, and measure its behavior as it under- 
goes a LZ transition. We purposely operate in a regime in 
which the decoherence time scale t^ is much shorter than 
the adiabatic passage, and in which T 3> A, where A is 
the tunneling amplitude. While evidence of LZ transi- 
tions in sup erconducting qubits has been reported before 
id I 111 . II4I , to our knowledge this is the ffist direct mea- 
surement of LZ transitions in a superconducting qubit in 
the high T and strong decoherence regime. 

In the original LZ problem, the system Hamiltonian is 



Hs^-iAa^ + ea,)/2, 



(1) 



with e = vt, where cr^^z are Pauli matrices and ly is the 
sweep rate for the energy bias. We take |0) and |1) to 
be eigenfunctions of CTz, denoting the "left" and "right" 
states in a double-well potential which can represent the 
two flux states in a superconducting flux qubit. If at i = 
—00 the system starts in state |0), then the probability 
of finding the system in same state at time t — -|-cxd is 
exactly given by [3|, |j| 



Plz = e--^'/2^ 



(2) 



If Hamiltonian ([T]) describes the dynamics of the two 
lowest energy states in a multi-qubit adiabatic quantum 
computer close to the energy anticrossing, then ([2]) is the 
probability of failing to reach the final ground state in 
the decoherence-free system. 

Now suppose that the qubit is coupled to an environ- 
ment. The total Hamiltonian H — Hs + Hs + Hiat, com- 
prises the system (jlj and environment Hb parts, and an 



interaction Hamiltonian 



(MRT) peak are related by [K 



-ffint = -QfTz/2, 



(3) 



W^ = 2Ter 



(6) 



that provides coupling between the qubit and an operator 
Q that acts on the environment. Here, we only consider 
longitudinal coupling to the environment which repre- 
sents flux noise affecting the flux bias in a flux qubit. 
We don't specify Hb explicitly, because if environmen- 
tal fluctuations obey Gaussian statistics, then all aver- 
ages can be expressed in terms of the spectral density 
S(uj) = /_!^di e^'^*(g(i)g(0)). Here, (...) denotes aver- 
aging over environmental degrees of freedom. Hamilto- 
nian ([3]) is what one expects for the effective interaction 
Hamiltonian for a large-scale AQC at the anticrossing, 
regardless of the type of coupling of individual qubits to 
the environment [9|. 

An immediate consequence of coupling to the environ- 
ment is that the relative phase between the two terms in 
the wave function that correspond to the two energy lev- 
els becomes uncertain after some time, an effect known 
as dephasing or decoherence. Due to energy-time uncer- 
tainty, the dephasing time r^p is related to the uncertainty 
in the energy eigenvalues or so called broadening W of 
the energy levels: 1/r^ = W\e\/^/e^ + Is?. liW,T <^ A, 
then the system will have a well-defined ground state sep- 
arated from the excited state by a well-defined gap and 
thermal transitions will be suppressed. One would then 
expect that holds even in the presence of noise, al- 
though A may be renormalized by high frequency modes 
of the environment 13|, Il4j . The important question is 



now what happens when VF, T ^ A so that the broad- 
ened ground and first excited states merge into each other 
and thermal transitions completely mix them up. Here, 
we answer this question both theoretically and experi- 
mentally. 

In the regime A ^ W , the dynamics of the system 
become incoherent. Using second order perturbation in 
A and assuming that the environmental is dominated by 
low frequency Gaussian noise, the incoherent tunneling 
rate from |0) to |1) is given by [16[ 



roi(e) 



A2 
8 W 



exp 



(«■ 



H'-/^SH,., = P/ 



dw S{uj) 
27r to 



(4) 
(5) 



with backward transition given by rio(e) = roi(— e). 
The transition rates therefore exhibit a Gaussian peak 
with a center shifted away from the resonance point 
e = 0. The width of the transition region, W, which 
is a measure of the environmentally induced broadening 
of the energy levels, is thus given by the r.m.s. value 
of the noise. In thermal equilibrium, the width W and 
the position e^ of such macroscopic resonant tunneling 



These predictions have been experimentally confirmed 
using superconducting flux qubits 



m 



Let us now return to the LZ problem. Suppose att^ti 
the system starts from |0) with probability Po{ti) — 1. 
In the incoherent tunneling regime (A <C W), the off- 
diagonal elements of the density matrix vanish quickly 
(within time scale t^ ~ 1/W^)- Thus to find Po{t), one 
needs to solve the equation 



Po 



-FoiPq -I- rio(i — -Po) 



(7) 



In the low temperature regime T ^ W (but can be ^ 
A), Eq. © requires that W <^ 2ep, thus separating the 
peaks of Foi and Fio such that roi(— e^) <C roi(ep). One 
can therefore neglect the second term in ([7]); because 
at points where Fio is peaked, 1 — Pq ~ 0, and at all 
other points, Fio ~ 0. The probability Pq will then be 
approximately given by 






-(7rA^/2iy)K 



/27rW 



exp 



i^-ep? 



2W^ 



de. 



lie,. 



-oo, this equation becomes 



1 + erf f ^V^ 



(8) 
(9) 



(10) 



If also 6/ ^- oo , then k = 1, yielding ^ , which is exactly 
the LZ transition probability in a completely coherent 
system, in agreement with previous studies [15| |. If the 
condition W ^ T does not hold, then one must keep all 
of the terms in ([7]) and calculate Po numerically. 

We have experimentally tested the above predictions 
by examining LZ transitions using a single qubit in a 
28 qubit chip designed for adiabatic quantum computa- 
tion. The sample was cooled down in a magnetically 
shielded dilution refrigerator with heavily filtered lines 
to a base temperature of about 10 mK. The qubits on 
the chip were compound Josephson junction (CJJ) RF- 
SQUID qubits as schematically shown in Fig. [T|d and de- 
scribed in Ref.|l7|. Each qubit consists of a main loop and 
CJJ loop subjected to external fiux biases ^x and ^^^J , re- 
spectively. The CJJ loop is interrupted by two nominally 
identical Josephson junctions connected in parallel. This 
device can be operated as a qubit for <I>^JJ S [0.5, l]<i>o 
and $j; ~ 0, where $o is the fiux quantum. The two op- 
positely circulating persistent current states correspond 
to the states |0) and |1). The bias energy is e = 2 |Jp| $3,, 
where the Ip is the magnitude of the persistent current. 
The parameter A is the amplitude of the flux tunneling 



between the two states, 
by $J^JJ. Maximum A 



Both \Ip\ and A are controlled 
'^ Wi, ~ 20 GHz, where Wr, is 




FIG. 1: Schematics of (a) the pulse sequence for the LZ mea- 
surements, (b) a CJJ RF-SQUID qubit, and (c) the qubit 
potential during the measurement. 
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FIG. 2: LZ probability as a function of inverse sweep rate for 
different values of A (<1>S^^ from -0.925$o to -0.91<I>o). The 
lines are linear fits to the data. 



the plasma frequency of the RF-SQUID, is obtained at 
$^jj = $0/2. For $^JJ « $0 one expects A ^ and 
the system becomes localized in |0) or |1). One can then 
read out the qubit by measuring the flux via an induc- 
tively coupled DC-SQUID (not shown). 

We isolated one of the qubits by tuning the coupling 
between that qubit and its neighboring qubits to zero, 
which allowed us to perform single qubit LZ measure- 
ment using the pulse sequence shown in Fig. [T^. The 
qubit is first initialized in one of the states |0) or |1) at 
a bias ^x = ^xA where the transition out of the state 
is very unlikely, and then the bias is linearly swept from 
^x,i to a final value ^x.f, at which point the qubit is mea- 
sured. A cartoon of the qubit potential during the pulse 



FIG. 3: (a) First MRT peaks in Fqi and Fio and their best 
fit with Eq. Q. (b) The experimental value of n as defined 
in ([SJ as a function of the final bias e/. The solid line shows 
the theoretical curve (|10|) using the parameters obtained from 
the best fit in (a). 



sequence is shown in Fig. [T]:. The probability Po{tf) 
of finding the qubit in the same state |0) as it started 
from was measured by repeating the above process 2048 
times for each value of the sweep rate ly. Figure [2] shows 
the probability Po{tf) in logarithmic scale as a function 
of 1/u. The result shows exponential dependence upon 
1/z/, in agreement with Eq. ([8|). 

Next we experimentally verify equation (jlOp . We first 
determine A, W and Cp by measuring Fqi and Fio, as 
described in Ref. |l7l Figure [3^ shows example plots of 
Fqi and Fio as a function of bias ^x for the above qubit 
at ^^JJ = —0.749 $0- The line-shape of the resonant peak 
fits very well with the Gaussian function (j4]), providing 
A, W, and Cp as fitting parameters. For the data shown 
in Fig. [3Ji, we found A = 0.082 ± 0.002 mK, W^ = 123 ± 2 
mK, and e^ = 354 ± 3 mK. Equation ([6]) then gives the 
effective temperature of the sample to be T = 21 ± 1 
mK. Notice that the condition T <^W ^ 2ep is approx- 
imately satisfied and therefore ^ should be sufficient to 
describe the LZ probability. 

The LZ probability as a function of flux bias was then 
measured for the same CJJ setting as in Fig. [3^. Fig- 
ure[3jD shows k = — \nPo/{TTA/2v) as a function of ^xj 
using the extracted A. The data starts from zero where 
^xj ~ ^x,i and shows a plateau at k « 1, in agree- 
ment with the theory. We have also plotted, on the same 
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FIG. 4; (a) Plot of A vs. flux bias ^l^^ applied to the com- 
pound junction. The lower curve is A measured with MRT 
and the upper curve is A measured with LZ. (b) Experi- 
mental (dots) and theoretical (line) ground state probability 
Pg = 1 — Plz, for flxed i/ = 0.05$o/^s as a function of A. 
The dashed line indicates the crossover between the incoher- 
ent and coherent regime, defined hy W ~ A. 



graph, the theoretical prediction of Eq. PH)) . using ep and 
W extracted from the MRT in Fig. [2^, and found very 
good agreement with the experiment with no extra fit- 
ting parameters. At larger biases, \^x.f\ > 2.5 m$o, the 
experimental k deviates from the theoretical curve due 
to tunneling to the first excited state in the target well. 

The measurements of the transition rates and the LZ 
probability allow us to extract A as a function ^^'J for a 
large range of A [l8| . Exponential dependence on $|^JJ is 
evident in Fig. 3^. Figurejljj plots the LZ probability, for 
a fixed value of ;/, as a function of A for a quite wide range 
of A (from 27 /kK to 1.25 K) together with the theoretical 
prediction [l9| . In the figure we have identified a line A = 
W which separates coherent tunneling from incoherent 
tunneling regime. Excellent agreement with theory is 
observed. 

We have reported on an experimental probe of the 
practically interesting regime for adiabatic quantum 
computation, where the energy gap A is much smaller 
than both the decoherence induced energy level broad- 
ening W and temperature T. The method used iso- 



lates a single qubit in a larger-scale adiabatic quan- 
tum computer, tunes its tunneling amplitude A into the 
A <C W, T hmit, and forces it to undergo a LZ tran- 
sition. We find that the transition probability for the 
qubit quantitatively agrees with the theoretical predic- 
tions. In particular, we demonstrate that in this large 
decoherence limit, the quantum mechanical behavior of 
this qubit is the same (except for possible renormaliza- 
tion of A) as that of a noise- free qubit, as long as the en- 
ergy bias sweep covers the entire region of broadening W. 
The close agreement between theory and experiment for 
a single qubit undergoing a LZ transition in the presence 
of noise supports the accuracy of our dynamical models, 
including both the noise model and the model of a sin- 
gle superconducting qubit that has been isolated from its 
surrounding qubits in an adiabatic quantum computer. 
Future experiments will test the behavior of multiple cou- 
pled qubits undergoing a LZ transition in the presence of 
noise. 
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